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We propose a way to control the Josephson energy of a single Josephson junction embedded in one- 
dimensional superconducting metamaterial: an inhomogeneous superconducting loop, made out of a 
superconducting nanowire or a chain of Josephson junctions. The Josephson energy is renormalized 
by the electromagnetic modes propagating along the loop. We study the behaviour of the modes 
as well as of their frequency spectrum when the capacitance and the inductance along the loop are 
spatially modulated. We show that, depending on the amplitude of the modulation, the renormalized 
Josephson energy is either larger or smaller than the one found for a homogeneous loop. Using typical 
experimental parameters for Josepshon junction chains and superconducting nanowires, we conclude 
that this mode-engineering can be achieved with currently available metamaterials. 

PACS numbers: 74.50.+r,74.40.-n,74.81.-g,74.78.Na 


I. INTRODUCTION 

It is well-known that low-dimensional superconductors such as Josephson junctions^, thin films?, and narrow wires? 
sustain plasma excitations: low-frequency collective oscillations of the superfluid charge density. This is possible as the 
low dimensionality renders Coulomb interactions less effective, thereby reducing the characteristic plasma frequency 
to values well below the superconducting gap. As a result, at low temperatures, the damping due to quasi-particles 
is negligible. In extended homogeneous systems, these oscillations acquire a propagating character with a dispersion 
relation that is linear for one-dimensional systems and of the square-root type for two dimensional systems?. 

Propagating plasma modes have been observed in a variety of systems, including wire networks?, Josephson junction 
arrays??, and high-T c superconductors??. They play a fundamental role in the collective behaviour of low-dimensional 
superconductors. For example, propagating plasma modes provide the quantum fluctuations responsible for the 
occurrence of quantum phase-slips that eventually drive the superfluid-insulator transition?? - —. This is relevant in 
view of recent experiments probing quantum phase-slips in nanowire a 14 ' 15 and in Josephson junction chain o 16 ! 17 . The 
propagating plasma modes are also at the origin of the damping of vortex excitations in Josephson junction arrays??. 

Using modern fabrication techniques, the parameters characterizing the properties of the collective plasma exci¬ 
tations in low-dimensional superconductors can be chosen from a relatively wide range. This is especially true for 
nanostructured superconducting circuits such as Josephson junction chains and arrays. SQUID-based chains for in¬ 
stance sustain propagating plasma modes with a phonon-like dispersion, whose group velocity can be tuned in situ 
with the help of the applied external flux?. This fact has already been exploited successfully in various experiments 
where Josephson junction chains are used as metamaterials providing dedicated electromagnetic environments. Ex¬ 
amples are the use of Josephson junction chains to obtain a tunable environment in order to study the Bloch band 
dynamics of a single Josephson junctiona l 18 ! 19 and to observe the dynamical Casimir effect—. 

In this paper, we show how the controlled introduction of inhomogeneities in a superconducting metamaterial can 
enhance its use as a dedicated environment. In general, a spatial modulation of the metamaterial’s parameters will 
affect the propagating modes as well as their frequency spectrum. An example is the frequency shift found for the 
modes of a disordered Josephson junction chain, together with the formation of localized modes??. This offers the 
possibility of mode engineering: positioning the mode frequencies and controlling the mode amplitudes locally along 
the metamaterial in order to optimize its frequency response as well as the way it couples to the system one wishes 
to study. 

Specifically, we consider a one-dimensional superconducting metamaterial, forming a loop closed by a small Joseph¬ 
son junction. The loop is threaded by a magnetic flux which induces a persistent current in the loop, see Fig. [1] 
The propagating plasma modes along the loop give rise to quantum fluctuations of the phase difference across the 
small junction, thereby renormalizing the junction’s Josephson energy, and hence the persistent current response??. 
Comparing a homogeneous loop with a spatially, periodically modulated one, we show that the renormalization can 
be made either significantly weaker or significantly stronger, depending on the sign of the modulation amplitude, see, 
e.g., Fig. [5l We find that this effect is mainly due to the effect of the modulation on the behaviour of modes close to 
the junction. 
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II. THE MODEL 

We consider the system depicted in Fig.[TJ It consists of a one-dimensional superconducting metamaterial of length 
L, closed by a small Josephson junction to form a loop. The metamaterial can be characterized by a space-dependent 
capacitance per unit length c(x) with respect to ground as well as by a space-dependent kinetic inductance per unit 
length l(x). The junction has a bare Josephson coupling energy Ej$ and charging energy Ec = (2e) 2 /2 C, where 
C is the junction capacitance. The loop is threaded by a magnetic flux <f>. For later use we define the reduced flux 
/<j> = 27 tJ>/$o, where <f>o = h/ 2e is the superconducting flux quantum. 



Figure 1: (Color online) Single Josephson junction with Josephson energy Ej } o and capacitance C, embedded in a loop of 
length L made out of one-dimensional superconducting metamaterial. The metamaterial can be either a thin superconducting 
wire whose parameters (such as cross-sectional area or distance to a nearby screening gate) are spatially modulated, or a chain 
of Josephson junctions (see inset) with spatially distributed capacitances Cg } n and inductances L c h,n- 


A. Hamiltonian 


The low-energy properties of the system can be described in terms of the superconducting phase along the loop, 
4>(x), and its conjugate momentum, the space-dependent superfluid density II(:r), such that [cj){x) 1 n(a/)] = iS(x-x'). 
The Hamiltonian H describing the system then reads 


H = 



e c {x)li 2 {x) + ei(x) 


dcj){x) 

dx 



[e c (L)U(L) - e c (0)n(0)] 2 
Ec 


Ej, o cos [()> (L) - (j) (0)]. 


(1) 


The first term on the right-hand side describes the one-dimensional, inhomogeneous superconducting metamaterial; 
it is a sum of two contributions. The first one, quadratic in the superffuid density, corresponds to the electrostatic 
energy, where e c (x ) = (2e) 2 /2 c(x). The second one, quadratic in the difference of the phase gradient and the 
magnetic flux, corresponds to the inductive energy associated with the supercurrent in the metamaterial. Here we 
define ei(x) = (<f> 0 /27r) 2 /2Z(:r). Throughout this paper we will assume the phase fluctuations in the metamaterial to 
be small, so that nonlinear phenomena such as phase-slips can be ignored. This is achieved by imposing the condition 
e c < ei- For later use, we also define the space-dependent plasma velocity u p i(:c) = \J\/l(x)c(x) and the zero-frequency 
dimensionless conductance g(x) = y/c(x)/l(x)Trh/ (2e) 2 of the inhomogeneous metamaterial. The condition e c < e; 
then translates into the condition g > 1. The quadratic Hamiltonian describing the metamaterial constitutes a low- 
energy description, valid for energies smaller than some cut-off energy E max . The resulting long-wavelength theory 
describes spatial variations of phase and density along the metamaterial on length scales larger than the cut-off length 
Lq — hv p i/ E max . 

The remaining terms in Eq. 0 describe the charging energy and the Josephson energy of the single Josephson 
junction, respectively. For simplicity, we will neglect the junctions’s capacitance C throughout this paper. This is 
possible as long as the capacitance of the metamaterial close to the junction is larger than C. The smallest part of 
the metamaterial involved in the plasma oscillation is given by the cut-off length Lq = hv p \/E msix . Then we can set 
C = 0 if Lqc C or, in other words, when gEc E m&x . 
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The parameters E m ax , c, and l can be related to the parameters describing the actual system realizing the meta¬ 
material. If it is made out of a narrow superconducting wire with cross-sectional area S = w 2 embedded in a 
dielectric material with dielectric constant e and sitting at a distance d from a screening gate, 1/c = 2ln(d/w)/e and 
l/l = 2 n s e 2 S/m, with e the electron charge, m the electron mass, and n s the density of the superconducting conden¬ 
sate. The quantities c and l become space-dependent if one modulates the cross-sectional area, S(x), and the distance 
to the screening gate, d(x). The cut-off energy E max is given by the superconducting gap A. The metamaterial can be 
considered one-dimensional as long as the cross-sectional dimension w is smaller than the penetration depth, w < Xl, 

where Xl = J mc u ght /(47rn s e 2 ) with cn g ht the speed of light. 

The metamaterial can also be realized using a Josephson junction chain consisting of junctions with Josephson 
inductance L c h = ($>q/2tt) 2 / Ej tC h, where Ej tC h is the Josephson energy, connecting islands of linear dimension a 
and with ground capacitance C g . Then we have c = C g /a and l = L c h/a. The condition g > 1 assures that the 
characteristic frequency u>o = yJl/L c hC g < Ej^h/h, such that nonlinearities associated with the Josephson effect can 
be ignored. The parameters describing the chain become space-dependent by modulating the size of the junctions 
and the islands along the chain, such that C g —> C g ^ n and L c h —> L c h, n , see Fig. [I] We assume the chain’s modes to 
have frequencies well below the superconducting gap A of the islands. We also ignored the capacitance C c h of the 
junctions forming the chain. This implies that we assume the mode frequencies to be smaller than the Josephson 
plasma frequency fojJ p ^ c h = \/&Ej }C hEc,ch- The cut-off energy E max is therefore determined by the smallest of these 
two energies A and huip^h- 

We summarize the parameter identification for Josephson junction chains and superconducting nanowires in table|U 



JJ-chains 

superconducting nanowire 

inverse capacitance per unit length 

a/ Cg 

2 In (d/w)/e 

inductance per unit length 

Lch/ & 

m/(2n s e z S) 

plasma velocity u p i 

&/ \/L c hCg 

ciight (w/X L ) \/hi (d/w)/ne 

dimensionless conductance g 

VC g /L ch nh/(2e )* 

{w/8X l ) (ftcught /e 2 )y / ne/\n(d/w) 


Table I: Metamaterial parameters in terms of the parameters of JJ-chains and superconducting nanowires. 


B. Classical phase configuration 


To gain some insight in the behaviour of the system, we start by determining the classical phase configuration 
4> c i(x) along the loop, i.e., the configuration for which the sum of the potential energy terms of Hamiltonian (JT]) is 
stationary. It satisfies the differential equation 


d_ 

dx 


ei(x) 


f dcj> ci(s) 
\ dx 



( 2 ) 


implying that ei{x)[d(j) c \(x)/dx — f®/L] = A where A is a constant. Thus, the classical phase configuration is given 
by 


r x dx 1 

4>d{x) = </>ci(0) + A + fax/L. (3) 

Jo ei[x ) 

The constant A is determined by minimizing the total potential energy. Introducing <f> 0 = <p c i(L) — </> c i(0) = /$ + A/Ei, 
where we define the total inverse inductive energy of the metamaterial 

-1 =/_*£_ (4) 

Ei J ei{x)' 1 ’ 

0 


the minimization condition reads 


L* 

-7-(00 - U) + sin 0o = 0. (5) 

Here L*/L = 2Ej/E jq: in other words, L* is the length for which the energy of the supercurrents in the metamaterial 
and the Josephson energy of the junction are of the same order. 
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Figure 2: (Color online) Persistent current I/I c as a function of flux for a classical loop. Curves from top to bottom 
correspond to L*/L = 4 (solid pink), 1 (dashed purple), and 0.4 (dotted blue), respectively. 


Calculating the total potential energy U for the configuration © using the condition ©. one finds that it depends 
periodically on the magnetic flux /$. It is straightforward to calculate the persistent current /(/$) = (2e/h)dU/df<£ 
induced in the loop by the external flux. It is given by /(/$) = I c sin</>o(/$), where I c = 2eEjfi/h and </>o(/$) is found 
by solving Eq. ©. In Fig. [2] we plot the expected current - flux relationship for various values of the ratio L*/L. 
For short loops L*/L 1 the response is purely sinusoidal. Indeed, according to ©, </>o ~ in this limit and the 
phase-difference across the junction is completely determined by the flux threading the loop. In the opposite limit 
L*/L < 1, the persistent current shows a saw-tooth like behaviour as a function of /$. In this limit, the Josephson 
junction pins the phase difference (f >o to values close to integer multiples of 27 t. When /$ is increased from — n to 7r, 
the phase difference <po ~ 0 will increase slowly, cj )o = L*/$/L, causing a linear dependence on /$, /(/$) = L*I c f$/L. 
When /$ reaches the value 7r, cf>o jumps from the value L*7r/L to the value 27r — 7 tL*/L, causing the persistent current 
to jump from 7 tL*I c /L to —7 tL*I c /L. This jump is again followed by a linear increase /(/$) = i*/ c /L(/$ — 27 t). A 
similar jump occurs for /$ = — 7r. 


C. Quantum fluctuations 


We will now include the effect of the electrostatic energy stored in the metamaterial. This energy causes the 
phase to fluctuate around its classical value <f> c \. Setting </>(x) = (f> c \{x ) + xi x )> we obtain the effective Hamiltonian 
H e ff = Hq + Hj that governs the behaviour of the quantum fluctuations x(x), where 


H n = 


f 

Jo 


dx 


e c (x)n 2 (x) + e ; (x) 


dx( x ) 

dx 


Hj = -Fj,o{[x<T)-x(0)]sin^o+cos[</>o + x(T)-x(0)]}. 


( 6 ) 

( 7 ) 


In the remainder of this paper, we will treat the Josephson junction as a perturbation, assuming the Josephson 
coupling energy Ej t o to be small. This means that we focus on relatively short loops, such that L < L*. We diagonalize 
the unperturbed effective Hamiltonian H$ with the help of the mode expansions 

X ( x ) = \J\ K + 5 ( 8 ) 

n— 0 V 

Umax / *- 

2 H y 2^fx)^ n ^ _ a "] ' ( 9 ) 

where the cut-off value n max is defined through Juo nin ^ = E max . This yields the correct commutation relation 
[x (x), n (x')] = id (x — x') when imposing the completeness of ip n (x), 

Y. Ipn {%) j’n jx') =d(x-x'). (10) 

n 

The unperturbed Hamiltonian Hq takes its diagonal form 

Umax 

Ho — ^ ^ fcbJn (^n^n H"~ 1/2) , 

n= 0 


( 11 ) 
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provided the functions ip n ( x ) satisfy the inhomogeneous wave equation 


sje c (x) — \ ei (x) — \J e c (x)ip n (x) \ = -h 2 u>l^ n (x) /4 


dx . 


together with the orthogonality relation 


( 12 ) 



Ipn ( X ) 1p m {x) 


dx = 


(13) 


Finally, since no current flows through the ends of the metamaterial in the absence of the junction, we impose the 
boundary condition 


dipn (x) 
dx 


= 0. 

x=0,L 


(14) 


For later use, it is convenient to introduce the auxiliary functions 4 f „(a;) = y/e c (x)ip n (x). They satisfy the inhomo¬ 
geneous wave equation 


:(x) 


dx 


ei (x) 


d^n(x) 

dx 


= -h 2 U>l^ n (x) /A. 


The functions 'F n (x) are subjected to the weighted orthogonality condition 


/ [e c (x)\ (x) {x)dx = S n>m . 

Jo 


(15) 


(16) 


We can express the fluctuating part of the phase x( x ) directly in terms of the functions \F n , 


XI 


n —0 


\/Kr^ n ^ +a ”] ■ 


(17) 


D. Effect of quantum fluctuations: suppression of the persistent current 


In this paper, we will analyze the effect of the quantum fluctuations on the persistent current response of the loop in 
the situation where the Josephson junction can be treated as a perturbation. This is achieved by computing the average 
of Josephson part Hj , Eq. 0, of the effective Hamiltonian H e g with respect to the fluctuations \, Eq- governed 
by the unperturbed Hamiltonian Hq, Eq. ED- The average of the first term on the right hand side of Eq. 0 vanishes. 
The average of the second term implies computing (cos[</>o+x(U)~x(0)]) = cos (cos[%(£)—x(0)])- Therefore, (Hj) = 
Ej cos (fro , where Ej = £j,o(cos[x(T)-x(0)])- Using the identity (exp{z[y(L) - x(0)]}) = exp{-([x(L) — x(0)] 2 >/2}, 
we obtain Ej = E Jfi exp(- “ ^n/ 2 ) wit h 


— 


tvjj. 


-['Fn (L) - \F„ (0)] v / ^MM- 


(18) 


Here we used the fact that (a^a„) = [exp(/3feu n ) — 1] _1 with /? = l/fc^T. This implies that the persistent current 
response of the loop will be qualitatively similar to but quantitatively different from the one found above in the 
classical limit: the bare Josephson coupling Ej t o energy should be replaced by its renormalized value Ej. 

The renormalized Josephson energy Ej has been computed in Ref. {HI for a homogeneous metamaterial. Setting 
c(x) = cq and l(x) = Iq, Eq. (fl5l) reduces to the homogeneous wave equation 


2 d 2 A/ n (x) 
Vpl dx 2 


u 2 n A> n (x), 


(19) 


where v 2 ^ = 1/IqCq. Its solutions are 


^ n (x) = y/2 e c /L cos (k n x), 


( 20 ) 
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Figure 3: (Color online) Ratio Ej/Ej t o, obtained from Eq. (1231) . as a function of n ma x for a Josephson junction embedded in a 
homogeneous metamaterial. Curves from bottom to top correspond to g = 1 (dashed-dotted green), 2 (dotted pink), 4 (dashed 
purple) and 10 (solid blue). 


with e 0 = (2e) 2 /2co, k n = mr/L and u> n = v p \k n . This solution satisfies the weighted orthogonality condition (fl()l) 
and is compatible with the boundary condition (1141) . Taking the zero-temperature limit, T = 0, we find 

= (21) 

where g = \/ca/loTrh/Ae 2 . Using the results (1A1I) and (1A4II from Appendix lAl we find 

IT-max 

^2 X 2 n = (2/<?)[Inn max + 7 + In2 + 0(l/n max )], (22) 

n= 0 

where 7 = 0.5772... is Euler’s constant. As a result 

Ej = Ej j o(l/ 2 e 7 n max ) 1/,ff . (23) 

We plot the ratio Ej/Ej.q as a function of n max for various values of g in Fig. [3] Below we will see how this 
result is modified for an inhomogeneous metamaterial. We will see in particular that, depending on the modulation, 
the renormalization of Ej^ down to Ej can be either weaker or stronger than the one found for a homogeneous 
metamaterial. 


III. MODULATING THE CAPACITANCE 


In this section, we consider a metamaterial for which the capacitance c(x) is modulated along the loop, yielding a 
modulation of the charging energy e c (x). Setting ei(x) = ei = (&q/2it) 2 /2lg constant, Eq. (TTSl) takes the form 


d 2 ^ n ( x ) U } 2 


dx 2 


^lO*) 


V I ; n (x) = 0, 


(24) 


where we used the fact that ei(x)e c (x) = Kv p \{x)/2. For simplicity, we restrict ourselves to modulations such that 

de c ( x) 


dx 


= 0 . 


x=0,L 


Then the boundary conditions for fkn (x) are d^ n (x) /dx\ x —o^L = 0. 

To be specific, we consider a periodic modulation of the capacitance, such that 

c (x) = cq(1 — t cos k m x), 


(25) 


(26) 


where Co is the average capacitance per unit length of the metamaterial, k m = rrm/L with m an integer, and t is the 
relative modulation amplitude such that |i| < 1. The modulation period is thus given by 2 L/m. It is easy to verify 
that E51) yields a modulation of e c (x) that is compatible with Eq. (l?5l) . 

With the modulation Eq. (E51) . Eq. (HM1) becomes 

(p \I/ (t;) 

— d ^ 2 J + E n [1 - t cos(k m x)] 'S/n (x) = 0, 


( 27 ) 
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where we defined E n = (tjj n /v v i) 2 with v 2 ^ = 1 /IqCq. Equation (1271) is similar to the Schrodinger equation. However, 
here the eigenvalue E n multiplies also the potential tcos(k m x). For this reason, even for small t. the potential term 
can become comparable to the kinetic term for sufficiently high energies. For t = 0, we obtain the homogeneous wave 
equation m discussed above. 


A. Perturbation theory with respect to t 


We wish to determine the eigenfunctions T^x) and the eigenvalues E n of Eq. m■ Since the modulation term is 
periodic in x, we expect the spectrum to consist of bands and gaps. As long as the relative amplitude t is small, t 1, 
we can use standard perturbative methods to obtain E n and Since the perturbation mixes the unperturbed mode 
cos (k n x) with the modes cos(fc n _ m x) and cos(fc n +mx), its dominant feature is to open a gap oc |t| in the spectrum 
E n at the degeneracy point k n = 7r?n/2L = k m /2. Other gaps exist, but they scale with higher powers of t and hence 
will be ignored in the following. We therefore set 

'ftn (^-) — A n COs(k n x) T B n COS (k n — m x') T Cri COs(fcy,_j_ m x), (28) 

and substituting this into Eq. (EZD. we find the set of equations 

0, (29) 

0, (30) 

0. (31) 

(32) 

'■'n+m 

The denominator of this expression is never small, as both n, m are positive, so C n = 0(t). Its substitution into 
Eq. (29) produces a correction to E n which is 0(t 2 ), so we neglect it in the following. At the same time, the difference 

kn-m ~ E n in Eq.(30) may become small at n « to/ 2, so degenerate perturbation theory for A n and B n should be 

used, which gives the characteristic equation for E n 

{kl - E n )(k 2 n _ m - E n ) - t 2 E 2 JA = 0. (33) 




From Eq. © we see that 


A n (fc 2 — E r j) H— —(B n + C n ) — 
t -^A n +(k 2 n _ m -E n )B n = 
-^-A n + (k 2 +m — E n ) C n = 


C n = 


tEj 2 


/A _ h 2 


This implies that 


C k l + k 2 _ m ) ± J(kl - k 2 _ m ) 2 + t 2 klkl_ m 

En.± = - 1 ^-> ( 34 ) 

where we ignored corrections of order t 2 (note that the square root is of the order t if the difference fc 2 — fc 2 _ m is 
small). The spectrum is given by E n _ for n < to/ 2 and by E n ^ + for n > m/2. For n = m/2, we find the expected 
gap |f|m 2 7r 2 /(4L 2 ). The gap region extends over a range Sn ~ mt around n = m/2. 

Equation (1321) works well in the gap region, however, it has problems at large n. Indeed, while the numerator ~ n 2 t, 
the denominator ~ mn. Thus, for n > m/t , Eq. (1321) no longer represents a small correction, so the perturbation 
theory breaks down. This breakdown is a consequence of the feature of Eq. m, discussed above: the weak modulation 
is multiplied by the eigenvalue E n , so it does not represent a small perturbation when E n oc n 2 is large. Thus, modes 
with large n should be treated differently, which will be done in the next subsection. For the moment, we restrict our 
consideration to the modes with n <C m/t, and proceed with the calculation of Ej{t). 

We next determine the coefficients A n , B n , and C n . In the vicinity of the degeneracy point such that n is within a 
distance Sn of m/2, the modulation mixes predominantly the modes cos {k n x) and cos (k n - m x), hence the coefficients 
A n and B n are of order unity, whereas C n is small, of order t. Then we have, up to a global normalization constant, 
A n ~ sin 9 n and B n ~ cos 6 n , where 


P _ 1,2 

tan 0 n = 2— - 

tEn 


( 35 ) 






Away from the degeneracy point, only A n remains of order unity, whereas B n decreases and becomes of order t , of 
the same order as the coefficient C n . Therefore, up to linear order in t, we have A n = D n sin 9 n , B n = D n cos9 n 
and C n = D n tE n sin9 n /[2(E n — where D n is a global normalization constant determined by imposing the 

weighted normalization condition, see Eq. m, 

L 

[ e c(x)]~ 1 'b^ l (x) dx = 1. (36) 


We find, to order t. D n = \j2ejL[\ + t sin(20„)/4]. 

With the same accuracy we obtain A„, using Eq. (fl8l) and setting T = 0, 


An — 


g\JL 2 E n V 




sin0„ {(-1)" - 1} + COS + ((-I)” m - l) + 


tE n /2 
E n — k 2 +m 


sin0 n {(—l)" +m — lj 


(37) 


For t = 0 we recover the result (l2lll for the homogenous metamaterial. For the modulated wire, the renormalized 
Josephson energy can be determined numerically using ED- summing A^, as long as ?r max <C m/t. Interestingly, as 
can be seen from Eq. ED, the corrections to A^ are of order t 2 for odd to, whereas they are of order t when to is 
even. This implies that for even to, the renormalization of Ej can be either stronger or weaker than the one found in 
the homogeneous case, depending on the sign of t. As we will show below, this parity effect is general and holds also 
when n max > m/t. 
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Figure 4: (Color online) Ratio Ej(t)/Ej(t = 0), computed from Eq. (1411) . as a function of n max . Curves from bottom to top 
correspond to g = 1, t = 0.1 (solid blue), g = 4, t = 0.1 (dotted pink), g = 4, t = —0.1 (dashed-dotted green), and g = 1, t = —0.1 
(dashed purple). 


Analytical results can be obtained in the limit t <C 1 /to, when the width of the gap region vanishes, Sn <C 1. We 
then approximate E n — k sin0„ ~ 1, and cos 9 n — (tn 2 /2)/[nr — (n — to) ]. Taking to even, we find a nonvanishing 
value A„ for odd values 21 + 1 of n only, hence 


A 


2 _ 
22+1 


1 

g (21 + 1) 


t (21 + l) 2 

(21 + l ) 2 - (21 + 1 — to ) 2 


t (21 + l) 2 

( 21 + 1) 2 - (21 + 1 + to ) 2 


(38) 


This is valid except for the point n = 21 + 1 = m/2 where the gap in the spectrum opens up. At the gap, according 
to Eq. (1371) . X ^/2 ~ 1 /to, which we neglect assuming to 1. We then evaluate 


^max 

A 22+l = (2/ff)[ 1 og4 + 7+ (1 + t/2)loglmax], (39) 

2=0,22+l^m/2 

where we dropped contributions of order 1/m and l/l ma x — 2/n max and assumed n max to. Result (1391) is obtained 
using results (IA2I) , (IA6I) and (IA8D from Appendix |Aj together with the identity 

21 + 1 21 + 1 1/1 1 
(21 + l ) 2 - (21 + 1 - to ) 2 + (21 + l ) 2 - (21 + 1 + to ) 2 ~ 2 V 21 + 1 - m + 21 + 1 + to 



As a result 


Ej(t) = E , j(l/4e 7 ) 1 ^(2/n max )^ 1+t / 2 ^ s . 


(41) 
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For t = 0we recover the result for a homogeneous metamaterial. In Fig.[I]we plot the ratio Ej(t)/Ej(t = 0), obtained 
from Eq. SB, of the renormalized Josephson energies found for the modulated and the homogeneous metamaterial 
as a function of n max for t = ±0.1 and for g = 1 and g = 4. Indeed, depending on the sign of t, the renormalization 
of Ej for a modulated metamaterial will be either stronger or weaker as compared to the homogeneous case. Since 
these results are obtained in the limit t -C 1/m, where the width of the gap region vanishes, Sn -C 1, we conclude that 
the effects of the modulation on the renormalization of the Josephson energy are due to the linear in t corrections to 
the modes, and not due to the opening of the gap in the spectrum. This is confirmed by the fact that the difference 
between the modulated and homogeneous metamaterials grows with n max : the larger the number of modes involved 
in the renormalization, the larger the effect of the modulation. 


B. The WKB approximation 


As we have seen above, the perturbation theory with respect to t breaks down for modes n > m/t. On the other 
hand, for modes with large n we can use a quasi-classical procedure to obtain the spectrum and the eigenfunctions 
of the wave equation (l27l) . In this subsection, we use the WKB approximation!! and solve Eq. (1271) for large energies 
E n . In this approximation, the mode (x) can be expressed as 


®n (x) 


\/Vn (X) 


M x ') dx ' + G n e-^o Pn (*')*’') , 


where p n is the quasi-classical momentum, 


(42) 


Pn (x) = V E n (1 - t cos k m x). 

The coefficients F n = G n in order to satisfy d^f n /dx\o = 0. Thus 


(43) 


(x) = 


\JVn (X) 


zF„ COS 


/ Pn{x')dx’. 

Jo 


(44) 


Imposing d^> n /dx\L = 0, we obtain the quantization condition f 0 dxp n (x) = mr. Thus the eigenvalues E n satisfy 


\[En = 


77 7T 2 t 2 

- - n7r -» k n [ 1 ± — ± 

2Ly/l^tE[2t/(t-l)\ 1 16 V >h 


where E(k) is the complete elliptic integral of the second kind, 

E (k) = f 2 y/l 
Jo 

and where we used the asymptotic expression^! 


— k sin 2 xdx. 




1 \ 2 k /1 ■ 3 \ 2 k 2 /1 • 3 • 5 \ 2 k 3 


1 2/1 V 2•4 


2-4-6 


y + ^(fc 4 ) 


(45) 


(46) 


(47) 


valid for k <C 1 , to obtain the approximate result for \JE n valid when |t| <C 1 . 

Next we use the weighted normalization condition Eq. (13611 to obtain the constant F n . Imposing 


4F 2 


z'V E n 


cos 2 < \/E n / \Jl — t cos k m x'dx ' > y/l — t cos k m xdx = 1, 


(48) 


defining s (x) = f* yjl — t cos k m x'dx' , and using the fact that s (0) = 0, s (L) = 2Ly/l — tE /vr, we hnally 

obtain 


4'n (x) = 


7re c 


cos y/E^,s{ x) 


Ly/l - tE (1 - tcoskmx) 4 


( 49 ) 


Note that both the energy eigenvalue Eq. (H5l) and the wave function Eq. (H51) are valid at all orders of t. 
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The validity condition of the WKB approximation used here reads dp/dx <C p 2 [23j. In our case, when |t| <C 1, 
this condition becomes n mt. On the other hand, as can be seen from Eq. (Tf5l) . the WKB eigenvalues E n do not 
reproduce the gap at n = m/2. We thus will use the WKB method only if 

n>m. (50) 

We see that, interestingly, there is a window m -C n <C m/t where both the perturbation theory and the WKB 
approach hold. As we will see below, this enables us to find the renormalized Josephson energy Ej in a broad range 
of parameters, combining both approximations. 

We finally turn to the computation of X n , Eq. m, using the WKB result (1491) . At zero temperature it reads 


fT 

(-If 1 

V gn 

_{1 - f (-l) m p (1 -t)*_ 


For t = 0, we recover the homogeneous result m ■ When m is an even number, we have 

\2 2 1-(—If 

" gn (1 - f)V2 ' 


(51) 


(52) 


This is of order t for small t. hence we expect that the renormalization of the Josephson energy depends on the sign 
of t in this case. When m is an odd number, starts from order t 2 . This parity effect is the same as the one found 
when using perturbation theory with respect to t. 


C. Renormalized Josephson energy 

We proceed and evaluate (cos[% ( L ) — x (0)]). We will combine the results obtained using perturbation theory (P) 
and WKB theory ( W ). We introduce the intermediate point n* = yjm ■ m/t = m/y/t within the interval m < n < m/t 
where both approaches are valid. Then 

(cos[y {L) — x (0)]) = e“^S XA » 1 »'e-iS=i( A »,p- A »^) 

= (cos[y (L) - x (0)})w ( cos bt ( L) - x (0)]) P _ W • (53) 

First we consider (cos[x (L) - x (0)]) w = exp[-(l/2) X l,w\ = (1/2e 7 n max ) 1 / 9 ' /T:ri , where we used JXl]) 

and (IA4I) . Compared to the homogeneous case, we see that g has been replaced by gy/l — t. Next we consider 

(cos[x (L) - X (0)])p_w = ex p[—(1/2) En=i ( X l,p - A U,1 v)]> where x l,p is S iven by Eq- (E3 and X 2 n W by Eq. (JH. 
In the range where the perturbation theory and the WKB approach are both valid, we use the latter with order t 
accuracy, hence 


^ +1 ’ w ~ g (21 + 1) 0 + l) ■ (54) 

The various sums can be computed numerically, some results are shown in Fig. El where we plot In Ej(t)/Ej(t = 0) 
as a function of m for n max = 100, g = 1 and t = ±0.1. First of all, we see that the modulation modifies the result 
for Ej significantly (by about 20 %) as compared to the homogeneous case. Also note the dependence on the sign 
of t. The effect is strongest for slow modulation; as m increases, the effect of the modulation weakens. This can be 
understood as the result of an effective averaging: a modulation with wave number m affects only those modes 
whose wave numbers n are larger than m. For modes with wave numbers n < to, the effect of the modulation of the 
capacitance along the metamaterial just averages out. 

Analytical results can be obtained in the limit t <C 1/m. Calculating P using Eq. ([38]) . we find 

A 2 i+i,p ~ x \i+i,w = (2i/s) {l/[2 (2 1 + 1) — to] + l/[2 (21 + 1) + m] — 1/(2 1 + 1)} . (55) 

Using (1A2I) . (IA6I) and (jA8|) . we obtain 


e 


5 E,n=l(^n,p-^n,w) = ( m g7)2j 


(56) 










Ej(t=±0.1) 

Ej(t=0) 
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Figure 5: (Color online) Ratio Ej{t) / Ej{t = 0) as a function of modulation wavenumber m for a metamaterial with n max = 100 
and g = 1 whose capacitance is periodically modulated with amplitude t = ±0.1. Both numerical results (red dots and blue 
triangles) and analytical results (green stars and purple diamonds) are shown. 


where we dropped again contributions of order 1/m and l// max — 2/n max and assumed ?z max m. As a result, 


Ej(t) 
Ej(t = 0) 


m 


2 n 


max 


t 

2 g 


(57) 


This result is also shown in Fig. [5] for <7 = 1 and t = ±0.1 (green stars and purple diamonds). We see that the 
approximate analytical result is fairly accurate for values of m < 30. This is more or less expected, as the condition 
t < 1/m breaks down when m > 10. However, the deviations between analytical and numerical results remain 
relatively small over the plotted range. This confirms the fact that the effect of the modulation is mainly due to the 
renormalization of the mode wave functions and not due to the opening of the gap in the mode spectrum. 


IV. DISCUSSION 

We verified that the results found above are not specific for a capacitance modulation. Indeed, introducing a 
modulation of the inductance l rather than of the capacitance c does not alter our conclusions. In Fig. [6] we show 
the results for the renormalization of the Josephson energy Ej for a metamaterial with c = c0 constant and 1 /l{x) = 
(l/Zo)(l — cos k m x). Details of the calculations are presented in Appendix iBl The results are very similar to the ones 
presented in Fig. 0 both qualitatively and quantitatively. Specifically, we find that the approximate analytical result 
for this case coincides exactly with the one found for the capacitance modulation, Eq. <EZD- The numerical result does 
not deviate strongly from the analytical result as long as to < 30. 


Ej(t=±0.1) 
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* *«»t .. 
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Figure 6: (Color online) Ratio Ej(t)/Ej(t = 0) as a function of modulation wavenumber m for a metamaterial with n max = 100 
and g = 1 whose inverse inductance is periodically modulated with amplitude t = ±0.1. Both numerical results (red dots and 
blue triangles) and analytical results (green stars and purple diamonds) are shown. 

We also investigated the case when the capacitance and the inductance are modulated simultaneously, such that 
c(x ) = co(l — tcoskmx) and l/l(x) = (1/Zo)(1 — tcosk m x). Such a modulation leaves the plasma velocity constant, 
ripi(a:) = f p i = y/l/loCQ. It corresponds, e.g., to a periodic modulation of the lateral size of the islands of a Josephson 















12 


junction chain, which affects the capacitance to ground C g and the Josephson coupling Ej }C h to the neighbouring 
islands in the same way, the latter being proportional to 1 /L c h. Details of the calculations are presented in AppendixlGl 
the results are shown in Fig. [7] The results are qualitatively similar to the ones presented in Figs. [5] and [6] Note 
however the quantitative difference: on a logarithmic scale, the effect of the modulation is twice stronger when both c 
and l are modulated. This is because both modulations contribute equally to the renormalization of the eigenfunctions 
T n . Indeed, comparing the approximate analytical result (|C15ll with Eq. (l57l) . we see that the modulation amplitude 
t is multiplied by a factor 2. The fact of having two contributions to the modulation also causes the difference between 
the approximate analytical result and the numerical result to grow faster with increasing m as compared to Figs. [5] 
and [TO 


Ej(t=±0.1) 

Ej(t=0) 


1.5 

1.0 




• f=0.1 (numerical) 

a f=-0.1 (numerical) 

* f=0.1 (analytical) 

♦ f=-0.1 (analytical) 


io 20 30 40 so m 


Figure 7: (Color online) Ratio Ej(t)/Ej(t = 0) as a function of modulation wavenumber m for a metamaterial with n max = 100 
and g = 1 whose capacitance and inverse inductance are modulated simultaneously with amplitude t = ±0.1. Both numerical 
results (red dots and blue triangles) and analytical results (green stars and purple diamonds) are shown. 


In order to see whether the effects discussed here are experimentally accessible, we give some typical values for the 
relevant parameters characterizing currently available metamaterials. Superconducting nanowire s^ 4 " — are character¬ 
ized by an inductance per unit length of about 0.002 H/m. The capacitance per unit length is about that of the 
vacuum, c ~ lOOpF/nr, but it can be strongly enhanced when putting the wire on dedicated substrates. For instance, 
using SrTi0?r with e = 10 4 as a substrate, one can achieve c ~ 10 6 pF/m. We thus have a plasma velocity in the 
range u p i ~ 10 4 - 10 6 m/s and a dimensionless conductance g ~ 1 - 100. Assuming the gap to be in the 10 Kelvin 
range, such wires sustain about 10 6 - 10 8 modes per meter. For a wire length of about 10/im, we have 10 - 1000 
modes. With these values, a 10 % spatial modulation of the nanowire parameters on a gm scale would realistically 
yield results similar to those shown in Figs. [5j EH and [TJ Josephson junction chains 7,19 may well have lengths N up to 
several thousands of junctions. Due to the relatively small capacitance to ground (typically, the ratio C g /C c h — 0.01), 
the characteristic plasma velocities are relatively high, such that the number of available propagating plasma modes 
nmax ~ N^/Cg/C c h typically is only about 10 % of the chain’s length. For the same reason, g ~ 1 is relatively small. 
Nevertheless, with these numbers a 10 % spatial modulation of the chain’s parameters over a few junctions would 
again realistically yield results similar to the ones obtained above. 

In conclusion, we studied the renormalization of the Josephson energy of a small Josephson junction, embedded 
in a loop made out of a one-dimensional metamaterial. We found that a slow modulation of the metamaterial’s 
parameters such as its capacitance or its inductance per unit length has a strong effect on the renormalized Josephson 
energy. Depending on the sign of the modulation amplitude, the renormalization can be either significantly stronger 
or significantly weaker than the one obtained for a homogeneous superconducting loop22. The modulation affects 
the electromagnetic modes propagating along the metamaterial as well as their frequency spectrum. However, inter¬ 
estingly, the effect of the modulation on the renormalization of the Josephson energy is mainly due to the modified 
behaviour of the mode spatial structure close to the Josephson junction, than due to the modified frequency spectrum. 
This example illustrates how mode engineering with a superconducting metamaterial can be used to affect its coupling 
to another superconducting device of interest in a controllable way. 
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Appendix A: Some useful sums 

In this Appendix we present results for the sums encountered in the main text. We start from the elementary sum^ 


— Inn max + 7 -\- 0 (\/ u max ), 

n 


(Al) 


where 7 = 0.5772... is Euler’s constant. When summing over odd values of n only, we havi 




E 21 T 1 =(1/2)lnZ max + 7/2 +In 2 + 0 ( 1/1 

max) • 


1=0 


From EB we see that a sum involving only the even values of n yields 


(A2) 


£27 = d/anii max + 7/2 + 0(1/1 max) • 


Using (IA2I) and (IA3I) . we obtain 


\n+l 


,l max / -1 \7 

^ ^-2 -= In 2 + 0 (l/n max ). 


n 


We next consider the sum for integer m > 0 


'I-max ^ 771 1 -j 71 max 771 ^ 

-= - V' - + V' - = ln(n max - m) - lnm + 0[l/(n max - m), 1/m], 

n. — m ' p ' p 


n=l,riy£m p—1 p—1 

where we used EB- We use this result to compute 


E 


i 


Z=l, 2 Z-l+m /2 

Similarly, for integer m > 0, 
1 


2 ( 2 Z-l)-m 


77 m ax “1“ 771 


—;- = -Y-+ Y - = ln(n max + rn)-lnm.+ 0 [l/(n max + rn),l/m], 

n A- m L ' 7i ^' 7i 


n + m p P 

71=1 p=l p= 1 


such that 


E 


^ 2 ( 2 1 - 1 ) + : 


= (l/4)[ln(n max + m/4) - ln(m/4)] + 0[l/(n max + m/4), 1/m]. 


(A3) 


(A4) 


(A5) 


= (l/4)[ln(n max - m/4) - ln(m/4)] + 0[l/(n max - m/4), 1/m]. (A 6 ) 


(A7) 


(AS) 


Appendix B: Modulating the kinetic inductance 

In this Appendix, we consider the case when the kinetic inductance l(x) is space-dependent and the capacitance 
c = Co is constant. Equation m then reduces to 

+ uj 2 n ^ n (x) = 0, (Bl) 

where v^x) = l/coZ(x). Now we consider a periodic modulation for the inductance such that 1 /l(x) = (l/Zo)(l — 
t cos k m x), in analogy with the capacitance modulation (1261) . As a result, r’pi(zr) = Upj(l —t cos k m : r), where = 1/loCo. 
We then can rewrite Eq. (IB1I) . such that 

(1 - t cos k m x ) d + tk ™ sin + En ^ n ( x ) = °> ( B2 ) 


d_ 

dx 




dx 
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where we defined E n = (ui n /v p i) 2 . 

We first use perturbation theory (P) and repeat the steps outlined in Sec. IIII Al for the capacitance modulation. The 
modes 'F„(:e) satisfy expansion (j28]), where the coefficients A n , B n and C n are now found from the set of equations 

( E n — fc^) A n + (t/2)k n k n — m B n + (t/2)k n k n+m C n = 0, (B3) 

( E n — B n + (t/2)k n k n - rn A n = 0, (B4) 

( E n — fc 2 _|_ m ) C n + (t/2)k n k nPm A n = 0, (B5) 

obtained using (IB2I) . From this we find that the eigenvalues E n are given again by Eq. (1541) : for the coefficients we 

find A n = ^2ejL sin 0 n , B n = yj2e c /L cos 9 n and C n = -t s /2e c /L sin 6 n k n k n+m /[2(E n - k 2 +m )]. Here, the angle 
6 n satisfies 


tan0 n = 2 


7.2 _ p 

tx 'n—m 

tk n k n — m 


(B6) 


As a result, the coefficients A n? p read 

A n.P = 


\ 5= [—ft-tt-D- - !(+»=».. -1} 


tk n k n +m/2 


sin#„ {(—l)" +m 


-1 


(B7) 


g\fL 2 E n 

When to is an odd number, A 2 P is even function with respect to t, whereas for even m the correction is of order t. 
In the limit t -C 1/m, we can set sinffy ~ 1 and cos0 n cs — (t/2)n(n — m)/[n 2 — (n — m) 2 ]. Then we obtain 


\ 2 

a 2/+1,P 


g (21 + 1) 


1 - 


t (21 + 1) (2Z + 1 — m) t(2l + l)(2l + l + m) 


(21 + l ) 2 - (21 + 1 - to ) 2 (21 + l ) 2 - (21 + 1 + to ) 2 


(B8) 


The validity condition of the perturbation theory is again n -C m/t. For larger values of n we can use the WKB 
method ( W ) to obtain the modes flfyfy). In this approximation, the modes are given by Eq. (1421) . where the quasi- 
classical momentum is now given by 


Pn (a;) = 


E n 


1 — t cos k m x 

Imposing the boundary conditions cW/ n /(ia:|o,L = 0, we see that F n = G n and find the quantization condition 




mr 2 y/l — t 


3t 2 

2LK[2t/(t - 1)] ' J kn ^ ~ 16 + ° 


where we defined the complete elliptic integral of the first kind 


K ( k ) = 


1 


\J l — k sm" x 


zdx. 


and used the asymptotic expression^ for small values of the argument, 


K(k) 


1 + 


1-3 

2^4 


1- 3-5 

2- 4-6 


k 3 + O (fc 4 ) 


(B9) 


(BIO) 


(Bll) 


(B12) 


to obtain an approximated value for E n when |f| <C 1. Defining, r (x) = [(( dx 1 /y/1 — t cos k m x' we finally obtain 

T , , _ / 7re c \/l - 1 cos y/E^r(x) 

^ n [X) — 


LK [2t/(t 1)] (l — t cos k m x ) 4 


(B13) 


imposing the weighted normalization condition (1361) . For even to and to order £, the parameter is nonvanishing for 
odd n only and given by 


X 2 

A 2l+l,W 


g (21 + 1) 


1+ 2>’ 


(B14) 
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which coincides with the result (15H) found for the case of a modulated capacitance. 

The WKB approximation can be used for large n tm. However, as in Sec. uni we use it for n^> m in combination 
with perturbation theory to compute the renormalized Josephson energy. Introducing again n* = m/y/t , we compute 
(cos[x(L) — x(0)]>, Eq. (1551) . numerically with the help of Eqs. (IB7I) and (IB 141) . An analytical evaluation is possible 
in the limit t -C 1 /m, with the help of (IB8I) . In this limit, 

2(2/ + 1) — m + 2(2/ + 1) + m_ ’ 
which coincides with the result (l55l) obtained in Sec. IIII Cl Hence the ratio E j(tME j(t = O') is again given bv Ea. (EZD- 


X 2 - X 2 
a 2Z+1,P a '21+1,W 


2 1 


Appendix C: Modulating both the capacitance and the kinetic inductance 


In this Appendix, we consider a simultaneous modulation of the capacitance and the inductance, such that c(x) = 
cq(1 — tcosk m x) and 1/Z(x) = (1//q)(1 — tcosk m x). Then, Eq. (fl5l) reads 


d 2l ^ n (x) (x) 

(1 — t cos k m x) —jt-1- tk m sin k m x — -1- E n (1 — t cos k m x) fUn (x) = 0, 


dx 2 


dx 


(Cl) 


where we defined E n = (ui n /v p i) 2 , with x 2 j = 1 /ZqCo. 

Using Derturbation theorv (P) as in Sec. IIII Al we start from the exDansion (1551). Using Eq. (ED, this yields a set 
of equations for the coefficients A n , B n , and C n , 


( E n k n ) A n + — ( k n k n — m E n ) B n + — (k n k n -\. m E n ) C n = 0, 

( E n k n — m ) Bn + “ ( knkn—m En) A n = 0, 

(A"n ^n+m) ^n + {J^n^n+m Bn) A n 0. 


Up to corrections of order t 2 , the eigenvalues are given by 


E n ,± = 


[kl + k 2 n _ m ) ± \J(kl - k 2 _ m ) 2 + t 2 (2k 2 n k 2 n _ m - k^k n - m - k n k^_ m ) 


(C2) 

(C3) 

(C4) 

(C5) 


This yields the gap 2|t|fc^, 2 for n = m/2. Up to linear order in f, we have A n = y/2e c /L[l + t sin(20„)/4] sin$ n , 
B n = ^/2e c /L[l + f sin(20„)/4] cos 9 n and C n = yj2e c jLt{E n - k n k n+m ) smd n /[2(E n - k 2 +m )\, where 


tan#„ = 2 


E n - kl 


n rv n—m 


t(E n h n kn-m) 


As a result, 


^ n.P — 


gV^E, 


= [1 + f sin(20„)/4] sin6> n {(-1)" - 1} + cos9 n |(-1)" m - l| 

En k n k n -\-m 


■sin 0 n (t/ 2)- 


R — k 2 
J^n rv n +m 


|(_l)n+m_ij 


For even values of in and in the limit (<l/m this yields 

4 


X 2 

a 2Z+1,P 


1 + t +t- {2l+1) 


5 (2I + 1) L 2 (2/ + 1) — m 2 (2/ + 1) + m_ 


(C6) 

(C7) 

(C8) 

(C9) 


The coefficients are nonzero for odd values of n = 2/ + 1 only. 

When n m/t, perturbation theory breaks down, and we resort again to the WKB approach. We obtain the 
eigenfunctions 


'kn (x) = 


12e c cos {y/E^x} 

L sj\ — tcoskmx 1 


(CIO) 
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which satisfy the weighted normalization condition (1361) . as well as the boundary conditions dA> n /dx\o,L = 0, provided 


The coefficients X n ,w then read 


\f~En = 


nir 

~L' 



(-If 1 

V gn 

yi-t(-i)"* vi-t 


For odd n = 21 + 1 and to order t we have 


A 


2 

21+1,W 


4 

g{2l + 1) 


(1 + 1 ). 


(Cll) 


(C12) 


(C13) 


Note that, as compared to Eqs. (15H) and (IB14I) . the amplitude t appears with a factor of 2 here. This reflects the fact 
that the capacitance and inductance modulation contribute equally to the renormalization of the wave functions 
As in Sec. IIII Cl we use the WKB approach for n > m in combination with perturbation theory to compute the 
renormalized Josephson energy. Introducing again n* = m/y/t, we compute (cos[x(L) — x(0)]}, Eq. (|53|) . numerically, 
now with the help of Eqs. (IC8I) and (IC12I) . An analytical evaluation is possible in the limit i < 1/m. In this limit, 
we use (IC9I) and obtain 


A 


21+1,P 


-A 


At 


21+1,W 


g [2(2/ + 1) — to 2(2/ + 1) + m ^ ^ 

Note that this is twice the result (1551) found in Sec. IIII Cl reflecting again the equal contributions of the capacitance 
and inductance modulation. The ratio Ej(t)/Ej(t = 0) is then given by 


Ej(t) _ / m \ o 
Ej(t = 0) \ 2n max J 

where the exponent has doubled as compared to Eq. 
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